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ABSTRACT 

Notwithstanding the advent of the Gamma-ray Large Area Telescope, theoretical models predict 
that a significant fraction of the cosmic 7-ray background (CGB), at the level of 20% of the currently 
measured value, will remain unresolved. The angular power spectrum of intensity fluctuations of the 
CGB contains information on its origin. We show that probing the latter from a few tens of arcmin to 
several degree scales, together with complementary GLAST observations of 7-ray emission from galaxy 
clusters and the blazars luminosity function, can discriminate between a background that originates 
from unresolved blazars or cosmic rays accelerated at structure formation shocks. 
Subject headings: diffuse radiation — large-scale structure of Universe — gamma-rays:theory 



1. INTRODUCTION 

The existence of diffuse 7-ray radiation was established 
in the early 70's at the dawn of 7-ra y astronomy by 
the SAS 2 mission ("Fichtel et alj ll975l ) and further as- 
sessed by COS B (Bennett 1990[). The determination 
of an extragalactic component (heretofore, CGB for cos- 
mic 7-ray background), which relies on a model for the 
foreground Galactic emission, became possible only later 
with the advent of the Compton Gamma Ray Obser- 
vatory. At energies above 30 MeV the EGRET experi- 
ment onboard the CGRO measured a CGB characterized 
by a specific intensity /i,(e)=5.4(e/keV)~^'^ ph cm~^ 
s~^ keV^^ sr~^ or an integral in tensity above 100 MeV , 
7=1.45 X 10^5 ph cm-2 g-i sr'^ (jSreekumar et al.lll998f ). 
A recent reassessment of this measurement changes the 
spectrum from a straight power law to a convex one and 
reduc es the integral inte nsity to 1.15 x 10~^ ph cm ~^ s~ ^ 
sr-i ([Strong et al.ll2004l . but see also Keshet et aL l2004f ). 

The origin of the CGB is still subject of debate. One 
naturally expects a population of unresolved blazars, 
a class of AGNs known to e mit up to the high- 
est energies, to contribute to it (fPadoyani et al.l Il993t 
iStecker fc SalamonlfTOOl iMiicke fc Poh]|l2000D . Ouanti- 
tative estimates, however, remain uncertain and various 
approa ches suggest a contribution o f order 25-50% of the 
CGB (jChiang fc Mukheried 119981 : iNarumoto fc Totanil 
l2006f) . Another candidate mechanism is inverse Comp- 
ton (IC) emission on CMB photons from cosmic-ray 
(CR) electrons accelerated at structure formation shocks 
both a round galaxy clusters (GC) and cosmic fila- 
ments (|Loeb fc Waxmanll2000D . This process can con- 
tribute up to a fraction 20% of the CGB, without vi- 
olating the existing EGRET up per limit s on t he 7-ray 
emission from individual GCs (jMiniatil |2002[) . Nor- 
mal and/or starburst galaxies have also been proposed 
as substantial 7 - ray e mitters (Pavl idou fc BrianI l2002t 
[Thompson et al.l l2007f ). Contributions from 7 r°-decay 
due t o CR protons in the core o f GCs (e.g., iMiniatil 
l2002t IColafrancesco fc Blasil Il998l ) and annihilation of 
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dark matter are also possible (e.g., lAndo fc Komatstil 
l2006t ). Note, however, that except for blazars none of 
these proposed sources has been detected in 7-rays as 
yet. 

The Large Area Telescope onboard GLAST** will ad- 
vance our understanding of the origin of the CGB, by 
improving our estimate of the Galactic emission, by re- 
solving part of the extragalactic component due known 
7-ray populations and, possibly, by detecting new types 
of 7-ray sources. According to current theoretical pre- 
dictions, a significant fraction (> 25%) of the CGB will 
remain unresolved by GLAST. In this Letter we study 
the angular intensity fluctuations of this residual compo- 
nent as a means to discriminate among viable models for 
its origin. We consider the case of blazars and CRs accel- 
erated at structure formation shocks. We use numerical 
simulations to model the relevant quantities involved in 
the calculations, particularly the correlation function of 
the sources. We predict a substantial difference in the 
level of fluctuations in the two scenarios above, testable 
with GLAST from a few tens of arcmin to several de- 
grees. This experiment, together with improved deter- 
minations of the blazar luminosity function and 7-ray 
emission from GCs by GLAST, will provide a discrimi- 
natory test for the unresolved CGB in the GLAST era. 
CGB fluctuati ons due to structure shocks were also com- 
puted in Wax man fc Loebl (|2000f ) analytically, and for 
an emission sufficient to account for the whole EGRET 
CGB. Our work is complementary to that of Ando et 
al. who, with different methods, studied the fluctua- 
tions produced by (a) sources similar to those considered 
here but resolved by GLAST (|2006a[ ). and (b) unresolved 
emission from blazars and annihilation in dark matter 
halos (i2006b'). Finally, CGB flu ctuations were also stud- 
ied by IZhang fc BeacomI (120041) for MeV emission from 
SN la. and bv ICuoco et al.l ()2006l ) for nearby emission 
above >100 GeV. Note that the success of any of these 
experiments hinges on the ability to assess and, if nec- 
essary, separate out fluctuations of Galactic origin. This 
requires a study the luminosity function and statistical 
spatial distribution of Galactic 7-ray sources similar to 
the one carried out below for extragalactic sources. 

Throughout we assume a ACDM cosmological 

* http://glast.gsfc.nasa.gov/ 
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model with flm—0.'3, flA=0.7, flb—0.04:, and h = 
i?o/100kms-iMpc"^=0.7 (|Spergel et al.l[2003l) . 

2. MODEL 

Given the integral intensity I{v), above photon en- 
ergy eo and along a direction v, the angular two 
point correlation function characterizing the fluctuation 
SI{v) = I{v) — about the solid-angle-average value 
{I)n, is the ensemble average C{9)—{6I{v)6I{u)), where 
0=arccos(u • v). When the correlation length of the 
sources is small compared to the spatial-scale over which 
they evolve and for small angular separations, <Cl ra- 
dian, we can write (|Limbeilll95a lPeebleslll993h 

where (j{z))v is the redshift dependent photon emis- 
sivity of unresolved sources in units phs^^cm^'^ and 
averaged over the comoving volume V, and C{9) is in 
units ph^ s~^ cm""* sr"^. The source correlation function 
is f(r), (11 is the luminosity distance and r{y,z)={y'^ + 
(i^O^y^^, with dc the comoving distance. For sources 
characterized by a rest frame spectral luminosity, L^{e), 
and a comoving luminosity function, C-y{L-y,z), with 
L^=vLi, computed at Eq in units erg s~^, we have 

{j{z))= dL^C^{L^,z)Lph{>eo[l + z]), (2) 

Jo 

where, Lph{> £o)= /J^ de e^^Ly{e), is the integrated pho- 
ton luminosity above a threshold Eq. In Eq. ([2]), Lumiz) 
is the luminosity of a source at redshift z with pho- 
ton luminosity Lph{> eo[l + z])—A'!rd'j^{z)Fiim, where 
Fum{> Eq) is the limiting integrated photon flux above 
Eo where a source is resolved. If the spectral emissiv- 
ity is a power law with index a, Li,[e) oc £~", then 
Liim{z) = AirdKa - 1)(1 + z)°'-'^eoFii,n- The angular 
power spectrum of the fluctuations is given by the trans- 
form, Cg=2n d{cos9)Pe{cos9)C {9) , where Pg is Legen- 
dre's polynomial of order £. Point sources also contribute 
Poisson noise t o the power spectrum and fo r unresolved 
sources this is (jTegmark fc Efstathioul[l996l ) 

1 dv/dz ['^-^'^^ , 

The average emissivity and correlation function of 
the CRs component are obtained through a numerical 
simulation of structure formation that follows the 
acceleration, spatial transport and energy losses of CR 
particl es. Th i s calc ulation uses the same technique 
as in iMiniatil (|2002f) . but a volume eight tim es as 
large. Shocks are identified ([Miniati et al.l l200Cf ) and 
a fraction of the particles crossing them is accelerated 
to a power-law distribution in momentum space with 
log-slope determined b y the shoc k Mach number in 
the test particle limit ([Belli 119781 ). The acceleration 
efficiency, ba sed on a variant of the thermal l eakag e 
prescription ([Kang fc Jones! Il995t iMiniati et all |2001| ). 
when expressed in terms of CR to shock ram pressure 
is always <40% and <1% for protons and electrons, 
respectively. The 7-ray emitting particles are affected 
by energy losses, dominated by synchrotron and IC 



emission for the CR electrons, and adiabatic losses, 
Coulomb and inelastic p-p collisions for CR protons. 
The latter process generates secondary which are 
also followed by the simulation as the primary electrons. 
The simulatio n is carried out w ith the cosmological code 
presented in Rvu et ahl (I1993D i n cora b inatio n with the 
CR code described in IMiniatil ([200X1 . l2002f) . We use 
a computational box of 100 Mpc on a side, with 
a grid of 1024^^ cells and 512'^ dark matter particles. 
Momentum space is divided into five log-bins. We 
compute the emissivity as a function of cell position and 
redshift, j{x,z), due to IC emission from CR electrons, 
and due to the combined emission from 7r°-decay and 
IC from e^. We simply refer to the latter as hadronic. 
The volume average is {j{z))=V~^ ^^^y z). 
The two point correlation function is then given 
by, £,{r,z)^{27T)~^ J Ps{k,z){s'mkr/kr)4Trk'^dk, with 
Ps{k, z) the power spectrum of 5(x, z)=j(x, z)/ {j{z)) — l. 

We compute the correlation function of blazars in the 
following way. We first use a high resolution N-body sim- 
ulation and adopt a set of physically and observationally 
motivated prescri ptions to populate colla psed halos with 
radio galaxie s (seelDi Matteo et al.|l200'i. for details). As 
illustrated in lDi Matteo et al . ( 2004) the model of radio 
galaxies accurately describes the observed angular corre- 
lation function as well as the radio luminosity function 
up to redshift z ~3. We rely on this model to com- 
pute the source correlation function. We use the estima- 
tor, e(r,z)=(7V''^[r,z]/7V;^^„,„„[r,z]) - 1, where 7Vi(r,z) 
is the sum over pairs whose distance, dij, falls within r 

and r+dr, and -^random i^ same quantity but 
expected for randomly distributed galaxies. We use this 
estimator instead of one in which the pairs are weighted 
with the product of the pair luminosity because for the 
angular scales relevant here it gives fully cons istent but 
much less noisy results ([Di Matteo et al.ll2004l ). 

Following previous work, we describe the 7-ray lu- 
minosity function of blazars, C^, with the luminosity 
function of AGNs at different wavelengths. We con- 
sider two models: the Pure Luminosity Evolution (PLE), 
based o n the radio luminosit y fun ction, originally pro- 
po sed inlStecker fc~ Salamon (|1996t ) and recently revised 
in iNarumoto fc Totani ()20()6l ) to improve the redshift 
distribution of EGRET blazars; and the Luminosity- 
Dependent Density Evolution (LDD E), based on the X - 
ray luminosity function of AGN (.Hasinge r et al.ll 2005f ). 
and studied in detail in INarumoto fc Tota ni (2006). We 
set Cj{Lj, z)=Ks{Ls/ L^)Cs{Ls, z) where s= R, X indi- 
cates radio and X-ray respectively, and assume a simple 
linear relation, Lj=10'^'' L^, between the 7-ray and the 
radio/X-ray luminosity of the sources. The relevant pa- 
rameters are (kr, (?r)=(0.173, 3.28) and (kx, gx) = (5.11 x 
10~^, 3.80). The assumed 7-ray spectrum is a power law, 
oc i^^", with a « 2.2. We use the luminosity function 
thus built to estimate the unresolved average emissivity 
from with Eq. Note that given the assumed linear re- 
lations between radio. X-ray and 7-ray luminosity of the 
sources, ^ does not depend on wavelength, as expected 
in unified AGN models in which the difference between 
blazars and radio galaxies is only ascribed to the jet ori- 
entation with respect to the line of sight. 

3. RESULTS AND DISCUSSION 
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In the following we consider the integral photon in- 
tensity above eo=100 MeV and adopt a flux sensitiv- 
ity for GLAST, Fiimi> eo)=2 x IQ-^ph s-^cm-"^. We 
provide a crude but robust lower limit for the unre- 
solved CR contribution by setting z^in ~ (^^o /c) {L^i^^ [> 
eo]/47rFH.m)i/2 Eq. (P), where i;r^(>eo) is the lumi- 
nosity of the brightest simulated GCs. This Zmin is the 
redshift beyond which the brightest and, therefore, all 
simulated GCs would be unresolved by GLAST even if 
they were point sources. We find zum ~ 0.05. 

The redshift evolution of the average emissivity for dif- 
ferent sources is shown in Fig. [T] With the chosen accel- 
eration efRciencies, IC (solid) and hadronic (long-dash) 
emission contribute about 20% and 7% of the CGB, re- 
spectively, most of which will likely remain unresolved. 
While the hadronic emission originates mostly in cluster 
cores, the IC emission is equal ly distributed in shocks 
around clusters and filaments (Miniati '2002) . This ex- 
plains the different curve normalizations, despite the fact 
that the cluster emission from both processes is below the 
EGRET upper limits. The unresolved blazars emissivity 
is also plotted for the LDDE (short-dash) and PLE (dot- 
long-dash) models, contributing a fraction about 30% 
and 20%, respectively, of the CGB. Note the marked dif- 
ference between the two models, with most of the con- 
tribution arising below and above z ~ 0.5, respectively. 
The correlation function at various redshifts is shown in 
Fig. [21 for the IC (top), hadronic (middle) and blazars 
(bottom) case. Note the large difference in amplitude 
for the correlation function of different potential CGB 
sources, especially at distances of several Mpc. Note also 
both the power-law shape and normalization change with 
redshift of the correlation function of blazars, reflecting 
the underlying correlation of the host galaxies. 

Fig. [3] shows the power spectrum of fluctuations of 
the integral photon intensity above 100 MeV. There 
SIe/{I)n = [£{e+l)Ce/4Tr]^/y{I)n, is the logarithmic 
contribution to the relative intensity variation by a mul- 
tipole i, corresponding to an angular scale 6 ~ 180° /£. 
The strongest signal is due to IC from structure for- 
mation shocks (shaded area). This is due to the much 
stronger correlation of shocks compared to blazars and 
their higher emissivity compared to hadronic processes 
(cf. Fig. [Hand [2D. At i < 100, or 6* ~ 1°.8, the fluc- 
tuations produced in the IC emission scenario are about 
10% of the average intensity and significantly above the 
Poisson noise from potential unresolved blazars, which 
is very similar in both the PLE and LDDE models 
(dot curve). Such fluctuations should be measurable by 
GLAST, given its angular resolution of 3-4° at lOOMeV. 
The fluctuations increase at smaller angular separations, 
£ > 100, still accessible by GLAST at higher energies, e.g. 
a few GeV, where the angular resolution is ~ 0.°5. Note 
that our results for the relative fluctuations are quali- 
tatively consistent with, a l beit a factor two lower than, 
those in lWaxman fc LoeE ()2000l ). Also, the fluctuations 
produced by hadronic emission (long-dash) are smaller 
than those due to IC emission in direct proportion to the 
ratio of the emissivities of the two processes. 

For blazars, in both LDDE and PLE scenarios the Pois- 
son noise dominates the signal for ^ > 50. Appreciable 
signal from spatial clustering of the sources is predicted 
only in the LDDE scenario (short-dash) at ^ < 50. At 



£ - 50, accessible by GLAST at 100 MeV, the inten- 
sity fluctuations are expected to be at the level of a few 
per cent, well below the signal from structure shocks. In 
the PLE model the angular fluctuations due to spatial 
clustering are much lower and always below the Pois- 
son noise. The reason is that in this model most of 
the unresolved emission is produced at high redshifts, 
where the amplitude of the correlation function decreases 
(cf Fig. [2]) and a correlated region appears projected on 
smaller angular scales in the sky (at least up to z <2 for 
the assumed cosmological model). This, however, does 
not affect the Poisson noise which only depends upon the 
number density and luminosity of the sources. 

Note that the CGB integral intensity is /(e) oc 
e~°'' , ai > 1. The spectrum of IC emission from struc- 
ture shocks is at least as flat as that (Miniati 200_2.), so 
that relative to the average intensity, the intensity and 
fluctuations contributed by this process should remain at 
least constant as a function of photon energy. In addi- 
tion, the integral sensitivity of GLAST up to a few GeV 
also scales as oc as is the integral spectrum of blazars. 
This implies that the number of resolved sources should 
not change appreciably as a function energy and that the 
Poisson noise from unresolved sources should also scale 
with photon energy as the background intensity. There- 
fore, the power spectrum predicted in Fig. [3| should be 
roughly independent of photon energy, between 100 MeV 
and a few GeV. This implies that by using information 
at different energies GLAST should be able to probe the 
power spectrum of angular fluctuations on a significant 
range of scales, from £ ^ a few tens at 100 MeV and up 
to a ^ 2± a few hundreds at 2 Gev. 

Various sources of uncertainty affect the results pre- 
sented in Fig. [3] Given an observed average CGB, {I)n, 
the intensity fiuctuations Sli predicted for the IC emis- 
sion from structure shocks are directly proportional to 
the assumed efficiency of CR acceleration. This parame- 
ter is highly uncertain and an efficiency lower by an order 
of magnitude would render the IC and LDDE model pre- 
dictions indistinguishable. However, in this case the IC 
model for the CGB would be ruled out, as it would pro- 
duce less than 2% of it. The important point is that 
if structure shocks contribute significantly to the CGB , 
GLA ST direct observations of nearby GCs (|Mimatill200"l 
l2003f ) would determine the efficiency parameter within 
a factor a few of the value assumed here. The pre- 
dicted large intensity fluctuations would then provide 
a signature of the (much larger) unresolved IC emis- 
sion from structure shocks. The consistency between ob- 
served 7-ray emission from GCs and CGB angular fluc- 
tuations provides a test for the IC origin of the CGB. At 
the adopted numerical resolution the structure of strong 
shocks, where the IC e mission is produc ed, should have 
numerically converg ed (|Rvu et al.l[2003f ). Nevertheless 
an uncertainty, AA^, in the determination of the shock 
Mach number, Ai, changes the log-slope q of the CR 

2 

distribution function by Aq = —^j^AAI. The shaded 
area in Fig. [3| includes the range of fluctuations obtained 
for a AM due to a change in the postshock pressure 
AP2/-F2±30%. It also includes variations obtained when 
changing Zmin from 0.05 to 0.01, which are negligible for 
£ >102 and lower the signal by - 30% at £ ~ 10. Fi- 
nally, we find that the occurrence of strong mergers adds 
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a sampling variance on the correlation functions, caus- 
ing a factor 1.5 uncertainty (reduceable with a larger 
simulation box) in the predicted fluctuations. The ambi- 
guity in the predicted blazars contribution to the level of 
CGB fluctuations at large scales is represented by the two 
curves for the LDDE and PLE models. However, this will 
be largely improved as GLAST which will discriminate 
between the two models based on their very distinct pre- 
dictions fo r the blazar luminosity fu nction and redshift 
evolution (Narumo to fc Totanill2006f) . 
We conclude that measuring the power spectrum of 



intensity fluctuations together with the faint end of the 
blazar luminosity function and 7-ray emission from GCs, 
should provide a valuable test of consistency for the sce- 
nario in which the CGB is produced by either IC emission 
at structure formation shocks or unresolved blazars. 
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Fig. 1. — Comoving volume-average emissivity as a function of redshift for IC emission from structure shocks (solid), Tr^-decay and 
associated e^ IC emission (long-dash), LDDE blazar model (short-dash), PLE blazar model (long-short-dash). 
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Fig. 2. — Spatial correlation function. Solid, dot, short-dash, long-dash, dot-short-dash and dot-long-dash indicate, respectively, 
redshift 2=0.01,0.1,0.5,1,2,3 for IC emission from structure shocks (top), 2=0,0.1,0.5,1,2,3 for hadronic emission (middle) and 
2=0.1, 0.4, 1, 1.6, 2, 3 (bottom) for blazars. 
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Fig. 3. — Power spectrum of angular correlation function for: 7-ray emission from IC at structure shocks (shaded area), hadronic processes 
in cluster cores (long-dash), LDDE blazar model (short-dash) and PLE blazar model (short-long-dash) and Poisson noise from unresolved 
blazars (dot). 



